Filling and wetting transitions of nematic liquid crystals on sinusoidal substrates 
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Close to sinusoidal substrates, simple fluids may undergo a filling transition, in which the fluid 
passes from a dry to a filled state, where the interface remains unbent but bound to the substrate. 
Increasing the surface field, the interface unbinds and a wetting transition occurs. We show that this 
double-transition sequence may be strongly modified in the case of ordered fluids, such as nematic 
liquid crystals. Depending on the preferred orientation of the nematic molecules at the structured 
substrate and at the isotropic-nematic interface, the filling transition may not exist, and the fluid 
passes directly from a dry to a complete-wet state, with the interface far from the substrate. More 
interestingly, in other situations, the complete wetting transition may be prevented, and the fluid 
passes from a dry to a filled state, and remains in this configuration, with the interface always 
attached to the substrate, even for very large surface fields. Both transitions are only observed for 
a same substrate in a narrow range of amplitudes. 

PACS numbers: 61.30.-v,61.30.Dk,61.30.Hn,61.30.Jf 



I. INTRODUCTION 

Wetting on smoothly structured substrates show a rich 
phenomenology. For simple fluids, the wetting behavior 
on structured substrates presents a variety of phenomena 
[lJ-Q j f° r which several physical laws were theoretically 
proposed, as the Wenzel |4| or the Cassie-Baxter laws [5|- 
These laws extend Young's law [(| for planar substrates, 
including new substrate geometrical parameters. If these 
laws are justified in some situations 0, they may also 
fail in general theoretical or experimental conditions Q. 
In addition to the usual wetting transition, simple fluids 
at bulk coexistence may present other surface transitions 
as the filling transition [T^]. From a thermodynam- 
ical point of view, the filling transition in simple fluids 
always precedes a complete wetting transition. However, 
the filling transition may not exist for shallow substrates 
@, but if it docs exist, it always precedes the wetting 
transition. 

The existence of long-range order in complex fluids 
may alter the scenario depicted above (Tlj . In a nematic 
liquid crystal, there is an orientational (continuous) long- 
range order. However, the presence of substrates which 
favor specific orientations frustrate the liquid crystal ten- 
dency to align along a given direction. As a consequence, 
clastic distortions emerge, altering the subtle balance be- 
tween the different free-energy contributions which lead 
to the surface transitions on structured substrates. 

In this paper, we show that the filling-wetting sequence 
which is observed in simple fluids on sinusoidal substrates 
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may be deeply modified when we have a nematic liquid 
crystal instead of a simple fluid. Depending on the phys- 
ical elastic parameters of the nematic, several scenarios 
may occur: we may only have a wetting transition; or 
a filling transition only, preventing the subsequent com- 
plete wetting transition; or the double sequence as in the 
simple fluid case. Both filling and wetting transitions are 
only observed for the same geometry in a narrow range 
of amplitudes. 

This article is organized as follows: in Section II. A. we 
briefly review the main simple analytical results about 
filling and wetting transitions of simple fluids. In Section 
II. B., we extend these results to nematic liquid crystal 
systems. In Section III, we use the Landau-de Genncs 
model of nematics to study their wetting behavior on a 
sinusoidal substrate. At the end, we draw some conclu- 
sions. 



II. MACROSCOPIC ANALYSIS 

A. Thermodynamic description of filling and 
wetting transitions for simple fluids 

In this section, we review the simple description of fill- 
ing and wetting transitions for simple fluids, following 
the ideas presented in [l(|. To that end, let us consider 
a structured substrate, like the one presented in Fig. [1] 
The substrate is translationally invariant along z, with 
L z being the total length along this direction, and pe- 
riodic along the x-axis with a period L. The substrate 
profile can be described by a simple even function g{x). 
To simplify the analysis, we consider a monotonically in- 
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FIG. 1. Schematic representation of (i) dry, (ii) filled and (iii) 
wet states of a simple fluid (1) on a sinusoidal substrate (w) 
of wave length L with vapor (v) above. 



creasing function for < x < L/2. 

We suppose the substrate favors the liquid phase and 
consider three possible interfacial states: i) dry, ii) filled 
and iii) wet (see Fig. [1]). If the bulk free energies of the 
vapor and liquid phases are the same (so they can be set 
equal to zero) , and if we neglect interactions between the 
interfaces, then the energies associated with each state 
will be simply the sum of surface energy terms. Let us 
denote S as the whole substrate area, A its projection on 
the x — z plane and T the excess interfacial free energy. 
Then, the free energy of the dry state (see Fig. QJi)) is 



T {€) = So,, 



(1) 



where a vw is the vapor-wall (substrate) surface tension. 

The free energy corresponding to the filled interfacial 
state (see Fig. QJii)) is 



2x, 



= S(x c )a lw + (S- S(x c ))a vw + I -^J Aa lv (2) 

where oi w is the liquid-wall (substrate) surface tension, 
(j[ v is the liquid-vapor surface tension, and S(x c ) is the 
substrate area in contact with the liquid phase, which can 
be obtained in terms of the abscissa x c of the contact line 
of the filled region with the substrate as: 



S(x c ) 



IX Vi + gldx 



g%dx 



(3) 



where g x is the derivative of g with respect to x. The 
value of x c can be obtained by minimisation of the free 
energy Eq. © with respect to it. Making use of Young's 
law, which relates the surface tensions mentioned above 
with the contact angle of a sessile drop of the liquid in a 
planar substrate, On, 



a tv cos tin 



(4) 



it can be shown that the minimum energy condition is 
fulfilled when x c satisfies (lfjj : 



= 2a lv (l - y/l + g%(x c ) cos6» w ) 

COS tin 



2a lv 1 



cos i/j 



(5) 



showing that the filling region must make contact with 
the substrate with an angle ip — 6 n . For smooth sub- 
strates, this solution is a local minimum only if g X x{%c) — 
d 2 g/dx 2 (x c ) < 0. 



Finally, the free energy for the wet state (see Fig. 
[T](iii)) is given by: 



J*** = Sa lw + Aa h 



(6) 



Different transitions between the different interfacial 
states can be observed. The transition from the dry state 
directly to the wet state is possible when the energies of 
these configurations are the same, that is, 

- = = S(a lw - a vw ) + Aa lv (7) 

Using Young's law Eq. (J4j), it is possible to deduce from 
Eq. ([7]) Wenzel's law for the dry-to-wet transition: 

— cos n = r cos 6^ = cos 6 r = 1 (8) 

where r = S /A > lis the roughness of the substrate. As 
ti r — > 0, a wetting transition may occur. However, this 
transition may be preempted by the presence of filled 
states. The filling or unbending transition between a dry 
and a filled state may occur if the free energies of these 
interfacial states are the same: 



p.«) _jrW = = S(x c )(a lw -a vw )+A [ ^ ) a lv (9) 



2x< 



which leads to the equation 



11° V 1 + 9 2 Jx 



2Xr. 



■ cos t 



= r c cos 9^ = 1 



(10) 



where r c > 1. If this transition occurs at a lower tem- 
perature than the wetting transition predicted by the 
Wenzel law (i.e. for larger values of On, as On usually 
decreases with temperature), then a filling transition oc- 
curs and Wenzel's law does not hold. Under these condi- 
tions, macroscopics dictates that the wetting transition 
between a filled and a wet interfacial state must occur 
under the same conditions as in the planar substrate (i.e. 
On = 0). However, this prediction can be changed as in- 
teractions between the substrate and the liquid-vapor in- 
terface are taken into account. So, if there is a first-order 
wetting transition for the planar substrate, the wetting 
transition is still first-order but shifted to larger values 
of On (smaller temperatures) [nil] . On the other hand, for 
planar continuous wetting its location and nature is un- 
changed @. 



B. Thermodynamic description of filling and 
wetting transitions for nematic fluids 

If instead of a simple fluid under saturation conditions, 
we have a liquid crystal at the nematic-isotropic phase co- 
existence, some new phenomena will inevitably arise. Let 
us consider that the substrate favors the nematic phase at 
a particular anchoring condition. From a thermodynam- 
ical point of view, if the bulk free energies of the phases 
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(iv) 




FIG. 2. Schematic representation of the nematic wet state 
on a rough substrate of typical I at isotropic-nematic phase 
coexistence. 



are the same, the total energy of the nematic interfacial 
state has a new positive extra term Td > 0, correspond- 
ing to the Frank elastic energy of the nematic distortions 
imposed by the anchoring at the substrate. Figj2] shows 
schematically an example of a nematic phase in contact 
with a substrate with homeotropic anchoring conditions 
(however, the theoretical arguments presented here are 
very general, and independent on the particular anchor- 
ing conditions involved). The total free energy of this 
interfacial state can be written as 



AtJi 



?d 



(11) 



The elastic deformation creates an effective long-range 
repulsion between the substrate and the interface: the 
closer the interface is to the substrate, the more con- 
strained the order will be, leading to higher energies. As 
the interface goes to infinity, the system will only have 
a characteristic length scale, that is imposed by the sub- 
strate (see Fig. [5]). If no defect arises, the elastic con- 
tribution may be simply estimated as Td = {KS jV)Tdi 
where K is the Frank elastic constant, and Td is a re- 
duced clastic contribution to the interfacial free energy, 
only dependent on the substrate geometry, but not on 
its scale I [l2j]. In some cases, topological defects nucle- 
ate close to or on the substrate. The presence of defects 
imposed by the substrate introduce a ln(//^) dependence 
on the reduced elastic contribution Td [l3[ . This depen- 
dence is easily justified. If the defects are induced by 
the substrate, they will place themselves at a distance 
of order of I (the characteristic length of the substrate) , 
usually close to its crests and/or troughs. On the other 
hand, the cut-off length is given by the size of the defects, 
which scales with the correlation length £. 

The dry-to-nematic wet transition may be derived from 
the generalized Wenzel law Q: 

j#v) _ j#) = o = S(a nw - a iw ) + Aa m + SjT d 



that can be rewritten as 



r cos( 



For large enough substrate lengths, I ^ K/<Ji n: the ef- 
fects of elastic terms are unimportant, and Wenzel's law 
is recovered. Usually, K and <n n are not independent, 
and their ratio scales with the correlation length £. In 
fact, using Landau-de Gennes model for nematic liquid 
crystals (see next section), it is possible to show that 
K/<Ji n ~ 10£ (at coexistence temperature). 

For the simple fluid case (cos 6^ = 1/r < 1), the dry- 
to-wet transition for the rough substrate occurs always 
before the wetting transition for the planar substrate, as 
the surface field or the temperature is increased. In con- 
trast, for the nematic case, the wetting transition for a 
rough substrate may occur either before (if the roughness 
is more important) or after the wetting transition for the 
planar substrate (if the elastic deformations are more im- 
portant). These deviations with respect to Wenzel's law 
are enhanced if defects are present, as the elastic term of 
the generalized Wenzel's law (Eq. (|12p) will decay slower 
with increasing I. The sawtoothed substrate exemplifies 
this effect [Hj. 

We now turn to the case of nematic filled states. Once 
again, it is useful to first consider configurations with- 
out defects, which do not appear if the substrate is shal- 
low enough. Given the fact that the nematic- isotropic 



interface favors a particular nematic anchoring [141 ]. if 
the interface is near the substrate, it will induce strong 
elastic deformations, because the nematic director is con- 
strained to follow the anchoring conditions on both the 
interface and the substrate. Thus the elastic energy of 
the deformations may be large enough to prevent any fill- 
ing transition, even if it would normally occur in the case 
of a simple fluid. It should also be pointed out that the 
interfacial free energies of the filled and wet states are 
normally very close to each other [To| . and a very tiny 
perturbation can exchange their relative stability. 

If the substrate imposes defects on the nematic ma- 
trix, a different picture emerges. The case of a sinusoidal 
grating is studied in detail in the next section. If the 
substrate is rough enough, and if the nematic-isotropic 
interface favors for example homogeneous planar anchor- 
ing, then the filled state (with the interface placed in be- 
tween the crests of the sinusoidal substrate) , may always 
have a smaller energy than the wet state for in the latter 
case, a periodic array of defects is necessarily created at 
the top of each crest. 



III. THE SINUSOIDAL SUBSTRATE. 
NUMERICAL RESULTS 

In this section, in order to substantiate our arguments, 
let us consider a nematic liquid crystal in contact with 
a sinusoidal grating, characterized by a wavevector q = 
2n/L, and an amplitude A (see also Fig. [1]): 



g(x) = A(l — cosqx) 



(13) 



rKTd 



(12) 



At the substrate, the nematic molecules preferentially 
align homeotropically, i.e. perpendicularly to the sub- 



strate. The system is translationally invariant along the 
out-of-plane axis z and periodic along the x axis. Fi- 
nally we impose that, far from the substrate, we have an 
isotropic ordering. 



A. The Landau-de Gennes model 

In the Landau-de Gennes (LdG) model, both isotropic 
and nematic phases can be locally represented by a trace- 
less, symmetric order-parameter tensor with components: 



3 11 

Qij = ^Slmrij - -Sij] + -B[klj 



(14) 



where n, are the Cartesian components of the director 
field n, S is the nematic order parameter which mea- 
sures the orientational ordering along the nematic di- 
rector, and B is the biaxiality parameter, which mea- 
sures the ordering of the molecules on the orientations 
perpendicular to n, characterized by the eigenvectors 1 
and m. In our model, we will only consider in-plane de- 
formations, although out-of-planc or twist deformations 
may also be important (a twist instability may occur 
for particular choices of parameters [Hj|). In this case, 
n = (cos#, sin (9, 0), and the tensor order parameter has 
only three independent components (namely Qn, Q22, 
and Q12). 

The LdG free energy may be written as 



•7"LdG = 



Wiik + <Pd)dV ■ 



^surf 



(l.S 



(15) 



where <^>buik is the bulk free energy density, <j> e \ is the 
clastic free energy density, and </> SU rf is the surface free 
energy, defined as (l6j : 

a Tr Q 2 - 6 Tr Q 3 + c[Tr Q 2 ] 2 (16) 

(17) 



fbulk 



= -^duQijdkQij + -^-djQijdkQ 



ik 



^surf 



;iuTr[Q • Q surf ] 



(18) 



where a depends linearly on the temperature, b and c are 
positive constants, and L\ and L2 arc positive param- 
eters related to the elastic constants. If wc rescale all 
the variables as follows (T^: Q = 6cQ/fe, the positions 
f = r/£, where the correlation length £ is defined as 
i 2 =_8c(3Li + 2L 2 )/b 2 , and F LdG = 24 2 c 3 F LdG /fb 4 , we 
get Ado = /y(0buik + 0d) dV + fg&suri ds, with rescaled 
free energy densities: 



0buik - \t Tr Q 2 - I Tr Q 3 + ^ [Tr Q 2 ] 2 



1 



3 + 2k 



rf = -~w)Tr[Q • Q S urf] 



(19) 
(20) 
(21) 




FIG. 3. Plot of the free energy densities (per projected unit 
area) of different branches of (meta)stable states, as a function 
of the anchoring strength w, for a sinusoidal substrate with 
L = 10£, k = 2 and Aq = 0.628, 1.288, 1.307 and 1.885. 
Numbers correspond to the typical configurations shown in 
Figffl See text for explanation. 



Here r = 24ac/b 2 is a dimensionless temperature, k = 
L2/L1 is an elastic dimensionless parameter (for stabil- 
ity reasons, the elastic parameter is restricted to the val- 
ues, k > —3/2) and w = Wwc/b 2 ^ is the dimension- 
less anchoring strength. Hereafter, we will consider these 
rescaled expressions, so we will drop the tilde notation. 

In our model, we will place ourselves at coexistence 
temperature, r = 1, for which the bulk free-energy den- 
sity has two minima, corresponding to 0buik = for 
rescaled scalar order parameters Si = (isotropic phase) 
and S*n = 1 (nematic phase) . It is important to note that 
the order parameter S in the coexisting nematic phase is 
rescaled, so its value in real units is 6/6c, which must be 
smaller than 1 (typically w 0.4). If the elastic parame- 
ter k is positive (negative), the nematic prefers to align 
parallel (perpendicular) to a possible nematic-isotropic 
interface. Finally, Q SU rf defines the favored tensor at the 
substrate. We favor a homeotropic alignment of the ne- 
matic by setting Q SU rf = (3i/®i/— 1)/2, with v the normal 
vector to the substrate, establishing a direct connection 
to previous papers [181 - 1201 ]. 



B. Numerical procedure 

For every set of model parameters, we numerically 
minimized the Landau-de Gennes free energy, using a 
conjugate-gradient method. The numerical discretiza- 
tion of the continuum problem was performed with a 
finite-element method combined with adaptive meshing 
in order to resolve the different length scales [21j . Due 
to the translational symmetry along the z axis and the 
periodicity on the x axis, we restricted the minimisation 
to a section of the system perpendicular to the z axis 
and with a width along the x axis equal to the period 
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L of the substrate. We used periodic boundary condi- 
tions on the lateral sides, and the fluid on the substrate, 
which corresponds to the bottom boundary of our unit 
cell, has not an imposed ordering or anchoring, although 
a nematic phase with homeotropic orientation is ener- 
getically favourable. Finally, we impose different fixed 
boundary conditions at the cell top (which we place at a 
large height H ) in order to explore the different intcrfa- 
cial states we may observe. The value of H was varied 
until we found convergence in the free energy. In order 
to find the dry (D) or the filling (F) interfacial states, 
we fix a bulk isotropic phase on top of the cell. Under 
these conditions, different states may be obtained when 
we vary the initial condition (for example, by considering 
initial conditions with different heights for the nematic 
filling region). The complete wet (W) interfacial states 
are obtained by imposing a fixed nematic phase at the 
cell top with a homogeneous nematic director along ei- 
ther the x axis (W") or the y axis (W^) (see Fig. 3] for 
illustrations) . The free energy for the wet states is calcu- 
lated by adding to the numerically obtained free energy 
the contribution of an isotropic-nematic interface paral- 
lel to the x — z plane, L(Ji n , calculated also numerically 
(and in some cases analytically) and corresponding to 
the most favoured anchoring conditions (homeotropic to 
the interface for re < 0, planar for k > 0). In this way 
we consider the presence of a nematic-isotropic interface 
infinitely away from the substrate in the wet state, al- 
lowing us to neglect the elastic deformations which may 
exist between the top cell and the nematic-isotropic inter- 
face. The true equilibrium state will be the state which 
gives the least free energy at the same thermodynamic 
conditions. 



C. Numerical results 

1. Case of re > 

As an example, in Fig. [3] we plot the free energy per 
projected unit area / = J- /A corresponding to the most 
relevant interfacial states for L = 10£, re = 2 (favoring a 
parallel anchoring at the nematic-isotropic interface) and 
different values of Aq = 0.628, 1.288, 1.307 and 1.885. 
For the shallowest case A = £ we plot the branches cor- 
responding to D states (which is the most stable for small 
w), W 1 - states and W". No F states are observed. As it 
is expected, the W" branch has always larger free energy 
than the others, since for shallow substrates the elastic 
deformations are stronger for this state than for the W 1 - 
state. On the other hand, the D and W" branches inter- 
sect at a first-order wetting transition for % = 0.2056, 
so W 1 - branch is the most stable for w > ww (typical 
configurations for both D and W 1 - states are shown at 
FigH(l) andH(2)). 

For rougher substrates, the scenario is completely dif- 
ferent (see Fig. [3]). For Aq — 1.885 four different (lo- 
cally) stable branches are observed: D states (again the 




FIG. 4. Typical configurations for a sinusoidal substrate with 
L = 10£: from left to right and top to bottom, (1) a D state 
(re = 2, Aq = 0.628 and w = 0.2), (2) a W r± state for shallow 
substrates (re = 2, Aq = 0.628 and w = 1.0), (3) a F" state 
(re = 2, Aq = 2.199 and w = 0.4), (4) a W 11 state for rough 
substrates (re = 2, Aq = 2.199 and w = 1.0), (5) a W 1 - for 
rough substrates (re = 2, Aq — 2.199 and w = 1.0), and (6) 
a F ± state (re = -1/2, Aq = 1.257 and w = 0.4). Numbers 
correspond to those shown in Fig. [3j (Color online) 



most stable for small w), W 1 - states, W" states and a 
filled state where the nematic is planar at the nematic- 
isotropic interface. This new state will be denoted by 
F" (a typical configuration is shown in Fig. HJ-(3)). For 
small w, the nematic-isotropic interface is almost paral- 
lel to the x — z plane. As w increases, the position of 
the nematic-isotropic interface for the F^ state is lifted 
and curved slightly, and it eventually pins at the crests 
of the substrate. Regarding the complete wet states, the 
free energy of the W 1 - branch is always larger than those 
of the W" states (unlike for shallow wedges). Typical 
configurations are shown in Fig. H(4) (W") and [1(5) 
(W 1 -). Now roughness induces larger deformations in 
W 1 - , leading to a higher free energy. In both cases topo- 
logical defects may nucleate either close to the substrate 
troughs (W 1 - states) or the crests (W^ states). We can 
see that the F^ branch has always a smaller energy than 
the wetting states, so there is a first-order filling transi- 
tion between a D and F^ state at wf = 0.1693. However, 
we do not observe a wetting transition, even for large w. 
This latter feature can be rationalised by the fact that 
the -F" state does not present topological defects as does 
the Wll state. 

The crossover from the scenarios described above oc- 
curs at intermediate values of Aq. Above some value of 
Aq the branch of -F" emerges as a metastable branch. At 
Aq fa 1.288 a triple point occurs, since the D, F^ and 
W 1 - branches intersect at the same value of w = 0.291 
(see Fig. [3]). Above this triple point, the F" states 
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FIG. 5. Global phase diagram for nematic adsorption on sinu- 
soidal substrates with k = 2. Black lines correspond the phase 
boundaries between the D, W ± and F" states for L = 10£ 
(solid lines), L = 20£ (dashed lines) and L = 40£ (dotted 
lines). Gray lines correspond to the transition from dry to 
filled state (dashed line) and from filled to wet state (solid 
line) in the case of isotropic fluid. Inset: zoom of the phase 
boundaries around the triple point for L = 10£. 



branch crosses both the D and W 1 - branches, so filling 
and wetting transitions are observed for the same geom- 
etry. However, the value of w for the wetting transition 
has a steep increase with Aq, so the wetting transition is 
no longer observable for Aq larger than 1.33. 

The global adsorption phase diagram for k = 2 is 
shown in Fig. [5j The domains for the different possi- 
ble stable phases (D, W 1 - or F") are divided by first- 
order transition lines which meet at a triple point. Note 
that the fH-W -1- transition line is almost vertical, so the 
range in which both filling and wetting transitions are 
observed is quite narrow. As the substrate wavelength 
L increases, the importance of the elastic term of the 
nematic free energy diminishes, and the D-W and D- 
F" transition lines shift towards lower values of to, ap- 
proaching the values obtained for simple fluids within the 
macroscopic theory outlined in Section II. However, due 
to the presence of elastic distortions, larger deviations 
than for simple fluids can be observed in all the cases, 

For the D-W- 1 wetting transition we can estimate 
this deviation through the modified Wcnzel law predic- 
tion, Eq. (1121) . The different terms in this expression 
can either be evaluated exactly or estimated numeri- 
cally. Firstly, the nematic order parameter for the LdG 
model close to a planar substrate in the presence of a 
bulk nematic or isotropic phase can be solved analyt- 
ically (l3| . The exact solutions of the corresponding 
Euler-Lagrange equations for the order parameter pro- 
files, d 2 S/dy 2 = 2S- 6S 2 + 4S 3 , are given by 



S±(y) 



5±(0) 



(22) 



S±(0) + (l-S ± (0))exp(W2</) 
where S+(y) and S—(y) are, respectively, the nematic 



order parameter profiles, as a function of the distance to 
the planar substrate y, corresponding to a bulk nematic 
or isotropic phase in bulk, and S± (0) are their values in 
contact with the substrate. The surface tensions between 
the planar substrate and the nematic or isotropic phase 
are given by 



/ [0ci(S+) + 0buik(5'+)] dy - wS + (0) 
= ^(1 + 25 + (0))(l - S + (0)) 2 - wS + (0) (23) 

poo 

a iw = / [0d(Sl) + <?Wk(£-)] dy - wS-(0) 
Jo 

= ^(3-2S-(0))S-(0) 2 ~wS-(0) (24) 
6 

The nematic order parameter values at contact with the 
substrate can be obtained from the boundary condition 
(dS±/dy)(y = 0) = —w as: 



s±(o) = 



1 



1= f(^ ±1 



(25) 



The planar contact angle 9 n is obtained from Young's 
law 



cos 9^ 



(26) 



where ct/tv ~ 0.178 for k = 2. 

The roughness r of a sinusoidal substrate, is given in 
terms of the amplitude A and the wavenumber q as: 



2E(-{qAf) 



(27) 



where E(x) is the complete elliptic integral of second 
kind. 

Finally, for the elastic term we have to resort to some 
numerical calculations. As a first estimate, we may use 
Berreman's approximation Fd ~ AK(Aq) 2 q/A [l^]. Fig. 
|5] shows the comparison between the numerical and the 
modified Wcnzel law predictions for the D-W 1 - transition 
line when Berreman's approach for the elastic free-energy 
contribution is used. They are in qualitative agreement: 
for small L the wetting transition is shifted towards larger 
values of w, and as L increases, the transition line moves 
down, approaching Wenzel's prediction, Eq. ((5J. A sim- 
ilar behaviour was observed for the wetting transition 
of the saw-toothed substrate @, EH- However, there 
are quantitative discrepancies even for the largest sys- 
tem considered. The explanation for these discrepancies 
is that Berreman's expression is only accurate for small 
Aq and wL very large (strong anchoring conditions), but 
it overestimates the exact result of Frank theory 22j. So, 
we have estimated the elastic energy obtained by using a 
Frank-Oseen (FO) model with a Rapini-Papoular surface 
contribution, where the interaction strength (which de- 
pends on to) is obtained from a contraincd minimisation 
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FIG. 6. Location of the D-W wetting transitions for n = 2 
and different values of L/£ = 10 (circles and solid lines), 20 
(squares and dashed lines), 40 (rhombi and dot-dashed lines), 
60 (upward triangles and dashed-double dot lines) and 80 
(downward triangles and double dashed-dot lines). Symbols 
correspond to the numerical results from full minimization 
of the LdG model (filled symbols), and the FO model (open 
symbols). Thick lines indicate the predictions from the modi- 
fied Wenzel law, Eq. (|12fl , using Berreman's approach for the 
elastic term. Thin lines only serve as guides for the eyes for 
the numerical results using the same style as the thick lines 
for same parameter L. For comparison, the prediction from 
usual Wenzel's law, Eq. ([8]), is shown (thick dotted line). 



of the LdG model in presence of a planar wall and a given 
nematic director orientation at a distance 1.5£ [l3j]. We 
used a similar adaptative-mesh finite-element method to 
minimise the FO model to that explained above for the 
LdG model. Our results show that, for Aq ~ 1 and 
assuming strong anchoring conditions, the elastic contri- 
bution is reduced, taking a value of around 75% of that 
obtained from Berreman's expression. However, if weak 
anchoring is considered, a further reduction is observed, 
even for L = 80£. This fact contrasts with our previous 
results in the sawtoothed substrate, in which the strong 
anchoring regime is reached for relatively smaller values 
of L. When the numerically obtained values of Td are 
used in Eq. (fT2|). the predictions for the wetting tran- 
sition from the generalized Wenzel law are in very good 
agreement with the full numerical values from the LdG 
model for L > 60£ (see Fig. [6]). 



2. Case of k = or n < 

Now we turn to the effect of n in the adsorption phase 
diagram. As it was already stated, for k — 2 the nematic 
fluid has a preferred planar anchoring at the interface. 
As k decreases the free-energy cost to anchor the ne- 
matic in a parallel or perpendicular orientation approach 
each other and become identical for k = 0. If we further 
decrease n, the system starts to prefer a perpendicular 



orientation at the interface. As a consequence, new in- 
terfacial states emerge, as the filled state F 1 - (see Fig. 
HJ-(6)). These new states appear at intermediate values 
of Aq for large w, and the NI interface shows a pro- 
nounced curvature following the substrate shape. The 
global adsorption phase diagrams for L = 10£ are shown 
in Fig. [7] Again, the continuous lines correspond to first- 
order transition lines which meet at triple points. For 
K = 0, the region of the phase diagram corresponding to 
F 1 - states is located between the regions corresponding 
to W 1 - and F^ states, with almost vertical boundaries 
between them. The D-F 1 - transition is weaker than the 
other transitions, so it is almost observed as a crossover 
for L = 10£. Finally, the VF" states free-energy branch 
approach to the free-energy branch. The dashed line 
corresponds to the (metastablc) D-W^ transition. For 
k = —1/2, the W" states become the stable states for 
large Aq and w. So, the phase diagram is reminiscent 
of that for k — 0, by swapping the F" states with the 
VF" states. It is interesting to note that it is possible to 
observe reentrance of F- 1 states as w is increased: the 
fluid undergoes a wetting transition from the F 1 - state 
to the VF" state, and for larger values of w there is a re- 
verse dewetting transition from the W" to the F 1 - state. 
However, this sequence of three transitions is limited to 
a very narrow range of values of Aq. 



IV. CONCLUSIONS 

In this article we have shown that the filling and wet- 
ting transition sequence of an isotropic fluid on a sinu- 
soidal substrate is deeply changed for nematic liquid crys- 
tals. Substrates favoring perpendicular anchoring in con- 
tact with nematic liquid crystals, that favor parallel an- 
choring at the nematic-isotropic interface, exhibit only a 
wetting (and not filling) transition for small Aq, while for 
large Aq only filling (and not wetting) transitions occur. 
In the latter case, the complete wetting transition is pre- 
vented by the pinning of the nematic-isotropic interface 
at the crests of the sinusoidal substrate, thus hindering 
the nucleation of topological defects that would other- 
wise be present in the wet configuration at large enough 
anchoring strengths. When perpendicular anchoring is 
favored at both the substrate and the interface, the filling 
and wetting transition lines change places. In this case 
the wetting transition may occur for larger Aq, whilst 
the filling transition occurs for intermediate Aq. Finally, 
in conditions such that the anchoring at the nematic- 
isotropic interface may be either planar or homcotropic, 
the phase diagram shows mixed features of the two pre- 
vious cases. The analysis of these situations is far more 
complex than in the previous cases, as any perturbation 
can change the delicate balance between the different 
contributions to the free energy 

In order to check the validity of the phenomenological 
generalized Wenzel's law, Eq. (fT2j) . we have compared its 
predictions with the numerical transition values for the 
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FIG. 7. Global phase diagram for L = 10£ and (left) k = ; 
(right) k — —1/2. The dashed line on the lefthandside figure 
is the D — W" metastable transition line. 

one-step D — W 1 - wetting transition. We see that they 
are in excellent agreement for L > 60£. This fact shows 
that elastic effects are important in the range of values 
of wavelength L considered in our work. In particular, 
they stabilize the D state with respect to the W 1 - state, 



as elastic deformations are a positive contribution to the 
surface free energy only present in the latter state. On 
the other hand, the location of the wetting transition 
approaches the value predicted by the typical Wenzel's 
law, Eq. ©, when L increases. 
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